Kinetics of a non-glauberian Ising model: global 
observables and exact results 

Sreedhar B. Dutta a , Malte Henkel 6 and Hyunggyu Park" 

"School of Physics, Korea Institute for Advanced Study, 
87 Hoegiro, Dongdaemun-gu, Seoul 130-722, South Korea 

b Groupe de Physique Statistiquc, Departement 1: Physique dc la Maticrc ct des 
Materiaux, Institut Jean Lamourjj] CNRS - Nancy Universite - UPV Metz, 
B.P. 70239, F - 54506 Vandceuvre les Nancy Cedex, France 
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behaviour after a quench to zero temperature is extracted. 
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1. Introduction 

Kinetic Ising models have been studied for nearly half a century, and have helped in 
addressing various fundamental issues in non-equilibrium physics. In spite of their 
simplicity, they have led to our understanding in various aspects, ranging from the 
relaxation mechanisms for equilibration to the collective behaviour in non-equilibrium 
systems. Unfortunately, even these simple models could only be solved exactly in very 
few cases, and one such case is the celebrated Glauber-Ising model [1]. 

The Glauber-Ising model is driven by spin-flip transitions and shows critical 
slowing-down near zero temperature with a dynamic exponent z = 2 in one 
dimension. Variational techniques and renormalisation-group methods suggest that 
spin-flip dynamics leads to z = 2, provided the energy-conserving flips are not excluded 
at any temperature [2j [3]. Kimball [I], and Deker and Haake [2] (KDH) studied a 
different kinetic Ising model, also driven by spin-flip transitions, but with rates which 
differ from those of the Glauber-Ising model, and find that z = 4 near zero temperature. 
Haake and Thol [5] further found a continuous set of single-spin-flip Ising models near 
zero temperature with z ranging from 2 to 4^ These models with z ^ 2 received a 
considerable amount of attention in the past, see e.g. [8] for reviews. 



| Laboratoire associe au CNRS UMR7198 

§ Their exact results are easily reproduced by heuristic arguments about the movement of domain 
walls [6] and have been extended to g-state Potts models [7]. 
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The equations of motion for total spin and three-spin in the KDH-Ising model are 
known to form a closed set of linear equations, and hence the global magnetisation could 
be evaluated exactly [H [2] . Since solvability for KDH dynamics is quite restricted, it has 
not been explored to the same extent as the Glauber-Ising model. Recently, two-time 
correlators and responses have been evaluated in the Glauber-Ising model to study its 
ageing properties P,[10l[TT] after a quench to temperature T = 0. It has been shown that 
the two-time space-time-dependent spin-spin correlator, C(t, s; r, r') = (a r (t)a r i (s)) , and 
the (linear) response function, R(t, s; r, r') = ^Ip^y , satisfy the scaling forms, if s 
and t — s > are sufficiently large (see e.g. [12] for a review) 

where a, b are ageing exponents, z is the dynamical exponent. Here spatial translation- 
invariance has been admitted. For y = t/s ^> 1, one further expects Fc,ii(0,y) ~ 
y-Ac,fl/2 ) where \c,r are the autocorrelation and autoresponse exponent, respectively. 
On the other hand, for correlations and responses with respect to the initial state one 
expects, for sufficiently large times t (see e.g. 



C(t,0;r,rO=r Ac/ ^c(^) , R(t, 0; r, r') = (^/'$ s {jqf) ( 2 ) 

Can one find similar non-equilibrium scaling forms in the KDH-Ising model ? And if 
so, what are the values of the exponents ? 

The layout of the paper is as follows. In the next section we define the model and 
derive the equations of motion for one- and two-point functions for spin and three-spin 
variables. In section [3l we derive exactly the susceptibilities and fluctuations of the 
global magnetisation and the global three-spin, and also discuss a dual description of 
the model as a reaction-diffusion process. We conclude in section HI In an appendix, 
details of the derivation of the equation of motion for global two-point correlators are 
presented. 



2. Kinetic Ising chains 

In this section, we define kinetic Ising models which evolve under single spin-flip 
dynamics, and explicitly find the equations of motion for expectation values and 
correlation functions for spin and three-spin observables. 



2.1. KDH dynamics 

We consider a kinetic Ising model on a chain A with energy 

H[(j] = - J 22 a nO- n +i - 2^ hn(Xn ' J > °> ( 3 ) 

where a spin configuration is denoted by a := {• • • , cr n _i, a n , a n+ \. ■ ■ •}, and a n = ±1 
is the Ising spin variable at lattice site n. The chain A C Z has L sites and we shall 
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take the thermodynamic limit L — > oo throughout. The system is assumed to be in 
contact with a heat-bath that induces only local spin-flip transitions, in other words, 
only single spins are flipped such that the transition rates depend just on the flipped 
spin a n and its nearest neighbours o n ±\. A configuration a, where the sign of the spin 
at site n is flipped is denoted by F n a := {■ • • , cr re _i, — & n , a n+ \. ■ • •}. We are interested 
in the following transition rates, whose form depends on the parameters 7, 5 

W(F n a\a) = a(l- |cx n (o" re _i + a n+1 ) + Sa^a^ ) 

- atanh(/3/i n ) (a n - ^ (cr n -i + a n+1 ) + <5<7 n _i<7 n cr n+ i) , (4) 

where (3 is the inverse temperature and a is a normalisation constant. For vanishing 
fields h n — 0, these are the most general local transition rates a 1— > F n o which are 
left-right (parity)-symmetric and also invariant under reversal under all spins [1]. The 
condition of detailed balance implies that 

7 = (1 + 5)tanh(2/3J). (5) 

The field- dependent terms are the most simple ones which are also compatible with 
detailed balance (for example, they are used in pQ [11], while in [9] a slight variant is 
studied). In what follows, we shall essentially restrict ourselves to the case 7 = 25 which 
we shall call KDH dynamics. Explicitly, when combined with ([5]), the value of 
7 _ tanh 2(3 J 

~ 2 ~ 2-tanh2/5J ' 1 ' 

and hence 8 — ► 1 when the temperature T = — > 0. As we shall see, this 
choice produces a closed set of dynamical equations for some global observables. This 
observation goes back to Kimball [1], and Deker and Haake [2]. The case with 5 = in 
equation (jlj) is the usual Glauber dynamics pQ. 

The master equation with only single spin-flip transitions is given by 

^P(a, t) = -J2 [W(F n a\a)P(a, t) - W(a\F n a)P(F n a, t)] , (7) 



and, following Glauber [T], can be rewritten as 
d_ 

Of 



Hence the dynamical equation for the iV-point function is given by 

d I N \ 

-Ql^rn ■ ■ ■ 0- nN ) t = -2 ^(T m • • • a nN ^ W { F nM ')l , ( 9 ) 

where n±, • ■ ■ , are N non-coinciding sites, and (■ • -) t '■— ^2 a ■ ■ ■ P{o~, t). 

Similarly, we can obtain equations of motion for two-time quantities from the 
conditional probability measure P(a,t\a, s) for a configuration to be a at time t 
conditioned to the configuration a at an earlier time s. The two-time correlation function 
of X : = 

o~m ■ ■ ■ o~ nN and Y :— a ni ■ ■ ■ o~ nM for t > s is defined as 

(X(t)Y(s)) := XYP(a, t\a, s) (10) 
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where Y := a ni ■ ■ ■ cr nM . The dynamic equation for P(a, t\a, s) is similar to equation (|S]) 
with P(a,t) replaced by P(a,t\a, s), and hence the equation of motion for two-time 
correlation functions is similar to equation (jHJ). 



2.2. Equations of motion 



Using equations (T4J) and ([9]), and re-scaling time such that a = 1/2, we get the following 
equation of motion for the expectation value of the spin: 



0_ 

di 



1 

(o~n) = - (o- n ) + 2 \o- n -x + °~n+i) - 8{q n ) 



+ (3h n (l - | (a n (o-^i + cr„+i)) + S(a n q n )j , 
to linear order in the field h n , where the three-spin variable q n is defined as 

q-n '■— o~ n -\a n a n+ i 
Similarly, the equation of motion for the average three-spin (q n ) is 



(11) 
(12) 



di 



{in) 



HQn) + l(o- n -i + o- n+ i) - 5(a n ) + 5 (A n ) 

p e » ' <- 1 



m=0,±l 

again to linear order in h n , and where 

7 



) + 5q n+m )^, (13) 



A 



25 



0~n-20~ n 0~n+l ~ O n -\0 n+ iG n+ 2 



+ 



_7 
25 



0~n~lO~ n O~n+2 ~ 0~ri-20~n-lO'n+l 



(14) 



In the absence of magnetic field h n , and under the assumption of translation invariance 
of the 3-point functions, (A n ) = if 7 = 25. Then ( TTTT) and (fTBl form a closed set of 
linear equations |U[2]. 

We now write down the equations for equal-time correlation functions when h n = 0. 
The evolution of spin-spin correlation function (o~ n o~ m ) for n 7^ m, obtained from 
equations (jlj) and (J9j), is governed by 



0_ 

di 



1 

(o~ m O- n ) = - 2(a m (T n ) + -(0- m {<7n-l + &n+l) + cr n (cr m _i + CT m+ i)) 



-8(o- m q n + a n q m ). (15) 

The correlation function (q n o~ m ), for |m — n| > 1, satisfies the dynamic equation, 

d 7 

^(om^n) = - 4(cr m g n ) + -(2cr m (cr n _i + cr n+ i) + q n (a m -i + cr m+ i)) 

- 5(a m a n + q m q n ) + 5(cr m A n ), (16) 

where A n is given in equation (1141) . The dynamics of the auto-correlations (cr n q n ) = 
(<7 n _i<7 n+ i) and that of the nearest-neighbour correlations (a n ±iq n ) = (a n a n ^\) are 
determined from equation (TT5]) . We further obtain the equation of motion for the 
correlation function (q n q m ), for \n — m\ > 2, to be 



9 I \ 

-^(QmQn) 



6{ 



iQn) + l(qm{<?n-l + 0"n+l) + <?n(o"m-l + 0"m+l)> 

l a n + q n a m ) + 5(q m A n + q n A m ). (17) 
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The dynamics of the nearest-neighbour correlations (q n <ln+i) = (c'n-i cr n+2) is given by 
equation ffT5]) . while that of the next-nearest one, (g n _ig n+ i) = (o" n _ 2 "n-i cr n+i cr n+2), can 
be obtained from equations (jlj) and Qj, and is given by 
d 

— (g n _ig rtf i) = - 4(q n - 1 q n+1 ) + 7(a-,« f icr„ +2 +c r n-20-„-i) + 25(cr n _ 2 ov 2 ) 

- 5(a n _ig ri+ i + g n _icr rtf i) + <5(g n _iA n+ i + g rM _iA n -i)- (18) 

Since A n appears in the equations (I13II15II16|I17II18I) . these equations are in general not 
closed and are therefore unsolvable. 



3. Global observables 



We now derive exact results for the behaviour of global magnetisation, response and 
correlation functions under KDH dynamics, when 7 = 25 is chosen. 



3.1. Magnetisation 

The global average of spin, M(t), and that of three-spin, T(t), are defined as 

M(t) ■= \ 5>»>* > T (t) := j (19) 

n n 

where L is the number of lattice sites. The equations of motion for these observables 
follow from equations fll l|13j) and read for h n — 

(20) 





Hence the global magnetisation M(t) and the global three-spin average T(t) are 

M(t) = ^ ((o+M(0) - 5T(0)) e" A -' - (a_M(0) - 5T(0)) e" A +') , (21) 
T(t) = ^ («_ (« + M(0) - 5T(0)) e- A - 4 - a + (e*_M(0) - 5T(0)) e"^*) , 



where 



A± = 2-5±A , A = (1 + 25 - 25 2 ) 1/2 , a± = 1 + 5 ± A. (22) 



In contrast to Glauber dynamics [I], the magnetisation can exhibit a non-monotonous 
behaviour due to the presence of two time-scales, see figured] It has a monotonous decay 
when the initial conditions are such that o^cT 1 < T(0)/M(0) < a+cT 1 . Otherwise, 
the global magnetisation first increases or decreases rapidly until a time t* and then 
decays more slowly with the finite time constant A_. Thus there is an initial increase 
in correlations with time before they begin to decay. We also see from figure [lb that in 
the zero-temperature limit (5 —* 1) the stationary value M(oo) no longer equals M(0) 
in general, in contrast to the case of Glauber dynamics. The crossover time is given by 



+> 

2A 



(a_M (0) - <5T(0))A„ 
(a+Af (0) - 6T(0))X. 



(23) 




Figure 1. Time-dependence of the global magnetisation M(t) in the ID Ising model 
with KDH dynamics for the initial values T(0)/M(0) = [0.50, 0.80, 0.90, 0.95, 0.99, 1.05] 
from top to bottom and for (a) S = 0.90476 . . . and (b) (5 = 1. The thick grey lines give 
the time-dependent global magnetisation of the Glauber-Ising model with the same 
values of (3 J. 



and depends on the ratio M(0)/T(0) and the temperature. This time-scale diverges 
either with fine-tuned initial conditions, a + M(0) — 5T(0) m 0, or by approaching low 
temperatures where 8 ~ 1. 

Now the equilibrium correlation-length £ at low temperatures behaves as 

r 1 = - lntanh(/5J) « 2e" 2/3J + -e" 6/3J . (24) 

3 

Using the equation (|24|) and the relation (jSJ) for 7 = 25, we get 

1 - r 2 + ^r 4 . (25) 

In this limit, for generic initial conditions which give rise to non- monotonous behaviour, 
the time-scale t* ~ 21n£ diverges logarithmically. Hence at low temperatures, the 
correlations gradually increase for a long transient time t*. After the crossover, for 
t ^> t* the relevant time-scale is the relaxation time associated with A_ and is given by 

r- = AI 1 « \i\ (26) 

and the dynamical exponent is [2] 

z = 4. (27) 

The reason for this exceptional value of the dynamical exponent is more transparent in 
a dual description, see below. 

Finally, we observe that if we choose 7 = —25, then the staggered quantities 

M s (t) ■= ~ E^ 1 )"^ , T s (t) := I (28) 

n n 
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process 


rate 


J 1 

dual 


process 


TIT — >TTT 


a(l + 7 + 5) 


AA 


^00 


TTT — >TIT 


a(l -7 + 5) 


00- 


— > 


TTi — -Tii 


a(l -5) 


0A - 


— > A0 


Til — -TTi 


a(l -5) 


A0 ■ 


— > 0A 



Table 1. Elementary spin-flip processes and their rates along with the dual 
description. 



satisfy the same system (1201) of equations of motion as those followed by the pair 
(M(t),T(t)) in the case of KDH dynamics with 7 = +25. Therefore, the conclusions 
reached for the ferromagnetic KDH model can be carried over to its anti-ferromagnetic 
analogue. 

3.2. Dual description: reaction- diffusion processes 

The nature of the relaxation dynamics along KDH line near zero temperature is more 
revealing in the particle-picture [13] of the Ising chain. In this dual description a particle 
(A) is associated with a kink (either 7 J, or j|) an d a vacancy (0) is associated with its 
absence (either TT or ID- Thus the spin- flip transitions are identified with the reaction- 
diffusion processes of the particles as shown in table [TJ A global reversal of the spins 
leaves the corresponding particle configuration unchanged, and hence only half of the 
spin-flip processes are listed in table [TJ 

Non- negative transition rates are found for the region ABC in the parameter space 
(5, 7) as shown in figure[2l The line AB is 7 = 1+5 and corresponds to zero temperature, 
while the 5-axis corresponds to infinite temperature. The line OB is the KDH line 
7 = 25, and as we move from O to B, the temperature drops from infinity to zero. The 
line BC is 5 = 1 where the diffusion process A0 « — > 0A is completely suppressed. Since 
the rate of diffusion D = a(l — 5), the particles or the kinks diffuse a distance £ in time 
r, where 

T = D = ^hry < 29 > 

As we approach the zero temperature the detailed-balance condition implies, as can 
be seen from equations (jSJ) and ff2M . that 7 « (1 + 5)(1 — £~ 2 /2). The parameter 
5 independently depends on temperature, and along KDH line near zero temperature 
5 ~ 1 — £ -2 . Now near the point (5,7) = (1,2) the most dominant process is the 
pair-annihilation A A — > 00, and for times larger than the corresponding transient 
time-scale the system reaches one of the many absorbing states of isolated particles. 
The remaining slow processes can induce diffusion and, from equation (1291) . we get the 
typical time needed to relax towards a state with domains of typical size £ as r m £ 4 [HI [5]. 
The nature of these dynamical classes can be qualitatively understood by comparing 
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Figure 2. The parameter space (J, 7) of spin-flip transitions is confined to the domain 
ABC; AB is T = line, OB is KDH line, OG is Glauber line, EC is 8 = 1, and OD 
is 7 = —2(5. 



the transition rates of energy-conserving flips A0 — ► 0v4 with energy-costing flips 
00 — ► For the KDH choice 7 = 25, the rates for energy-conserving-flips and 

for energy-costing-flips are the same ~ £ -2 , hence a large number of absorbing states 
are accessed, unlike the other near-equilibrium kinetic Ising classes with z = 2. The 
number of absorbing states for large L is of the order ((1 + v / 5)/2) i [T3j, while there 
are only two zero-temperature equilibrium states. 

A similar analysis can be made in anti-ferromagnetic Ising chains for the choice 
7 = —25. Along the line OD in figure [2] the staggered magnetisation can be exactly 
evaluated. The relaxation time r s for such a configuration at low temperatures can 
be easily estimated in the dual picture, and in ferromagnetic models r s ~ £ 2 , while in 
ant i- ferromagnetic models r s m £ 4 . 

3.3. Global responses and correlations 

We now evaluate the susceptibilities and fluctuations of total spin and three-spin. The 
responses of the magnetisation and of total three-spin to a uniform time-dependent 
magnetic field h(s) are defined as 

R{t,s) := — — ~FlTY ' Q ) ■= 777 X TT7T ■ 30 
PL ^ 6h(s) h=0 0L ^ 5h(s) h=0 

The fluctuations of the global spin and three-spin are described by 

Cm{t) := T^^,^n{n+m(t) > C 9f {t) := — ^ C%{ n (t), (31) 
n m,n 

where g and / are either a or q, and the equal-time correlation functions are given by 
C% n {t) := (a m a n ) t , C% n (t) := {q m a n ) t , C* B (f) := (g m g„) t . (32) 
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3.3.1. Response functions: The equations of motion for R(t, s) and Q(t, s) are obtained 
from equations fTTTj) and fl30l . They read 

^R(t, s) = (25 - l)R(t, s) - 6Q(t, s) (33) 

for t > s, while for t — s it is 

R(s, s) = 1 - 2£C? ff (s) + 5C™{s). (34) 

Similarly, we obtain the equation of motion for Q(t, s) from equations ([TBI and f[3"Uj) and 
have for t > s 

^Q(t,s) = 35R(t,s)-SQ(t,s), (35) 
and for t = s the initial condition 

Q(s, s) = 5 + 2(l-5)C™(s) + (l-25)C™(s) + 25C™(s) - 25Cf(s). (36) 
The equations fl33| 1351) are similar to (1201) . and we find 

R(t, 8 ) = ± (A4s)e- X -^ - A + (s)e- x ^) , (37) 

Q(t, s) = ^ (a_A_( S )e- A -(*-^ - a + A + ( S )e- A +(*- s )) , (38) 

where v4 T (s) = a-ti2(s, s) — 5Q(s, s) and a± are given in equation The initial equal- 
time responses R(s,s) and Q(s,s) depend on the functions C^(s), but for the chosen 
KDH dynamics, closed equations for them are unknown. In any case, the responses 
R(t, s) and Q(t, s) should decay, at large times, as exp(— (t — s)/r_). 

Aspects of non-equilibrium relaxation can be explicitly studied for responses with 
respect to the initial state, since closed expressions can be given for R(t, 0) and Q(t, 0). 
We study three examples: 

(i) Consider fully ordered initial states • ■ • TTTT ' ' ' an d ■ • ■ II! J, ■ • • at the critical 
point 5—1. Then both R(0, 0) = Q(0, 0) = and these global responses 
R(t, 0) = Q(t, 0) = of the system will not react to a small perturbing external 
field. Hence there is no perceptible equilibrium relaxation at criticality via global 
observables for ID KDH dynamics. 

(ii) In order to study non-equilibrium relaxation, consider a fully disordered initial state 
and quench at time t = the control parameter to the value 5. Then i?(0, 0) = 1 
and Q(0,0) = 5. For the limit 5 — > 1, we find R(t,0) = 1. This indeed describes 
a non-equilibrium relaxation, as may be seen by considering an initial state at 
thermal equilibrium with temperature Ti n i > 0. Standard techniques [15] may be 
used to calculate R(0, 0) and Q(0, 0) and we find that for 5=1 

R(t,0) = (1 -tanhn) 3 + tanhr/(l -tanhn) 2 e -2 * (39) 

where we have set rj := J/T^. The stationary value of the response function 
is distinct from the equilibrium value, which would only be reached if the limit 
rj — ► oo were taken. 
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(iii) In order to appreciate better the role of the initial state in non-equilibrium 
relaxation, consider the initial ensemble made from the states ■ • • TTITTITTI ' ' ' 
and • • • HT4ITIIT ' ' ' (with equal probability) and couple the system at the 
initial time to a heat bath such that the control parameter has the value 5. Then 
#(0,0) = 1 + 5/3 and Q(0,0) = -1 + 55. We find in the limit 5 -> 1 that 
R(t,0) = (4/3) exp(— 2t) relaxes exponentially fast towards its equilibrium value. 

For an interpretation in terms of dynamical scaling, we recall that the global 
responses calculated here are actually the Fourier transforms, viz. R(t, s) = 
J R dr e~ iqr R(t, s; r, 0) | _ Q of the space-time responses defined in the introduction, at 
vanishing momentum q = 0. Hence from (T5]) we expect R(t, 0) = t ( - 1 ~ XR ^ z J R du &r(u). 
For t sufficiently large, this may be compared in particular with equation ( |39i) . and this 
scaling interpretation suggests that 

\r = 1. (40) 

Provided that a + R(0, 0) — 5Q(0, 0) ^ 0, this result is generic. The independence of the 
initial state confirms the expected universality of the autoresponse exponent \r. 

3.3.2. Correlation functions: The equations of motion for the global correlation func- 
tions Cp \t) and C 9 f (i) are explicitly obtained in the appendix, see eqs. flA8IIA10|A12j) . 



While in general a closed system of equations of motion cannot be found, it turns out 
if the initial conditions C 9 ^(0) are of order 0(1), then in the limit L — > oo the global 
correlators C 9 '(t) satisfy the following closed set of linear equations 

C aa {t) \ 

-5 C q °(t) . (41) 
-6 / \ C"(t) i 




C qq (t) J 

Upon resolving, we find the global spin-spin fluctuations to be 

C aa (t) = B_e -2A -* + B e~ {x - +x+)t + B + e~ 2X+t (42) 
where A± are given in (I2"2"j) and 

B T := (a 2 ± C™(0) - 25a±C 9CT (0) + 5 2 C qq (0)) , (43) 

B ° ■= ~2^2 («+«-C CT<7 (0) - 25(1 + S)Cv{0) - 5 2 C qq (0)) . (44) 
The explicit expressions for the other two correlation functions are 

C q °{t) = \ (a_5_e- 2A -* + (1 + 5)B e^ x - +x+)t + a + B + e~ 2X+t ) , (45) 
5 

C qq (t) = 3 ( —B_e~ 2X - t + 5 e- (A - +A+)t + ^±B + e~ 2X A . (46) 
\ot+ a- J 

The global correlations at large times in general relax as exp(— 2t/r_). Since the leading 
correction, coming form the non-global terms, to these correlations will be of order 
0(t/L), the solution is valid for any time t < r_ in the large-L limit only for (B ,B±) 
of order 0(1). 
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The two-time correlation functions of a n and q n variables for t > s + is denoted by 
C% m {t,s) + := Y^9nLP(<r,t\v,s), (47) 

a, a 

where f n and g n are either a n or q n . These quantities have similar equations of motion 
as those of one-point functions of a n and q n . For any time t the correlation functions 
are defined as 

f C# m M + for t> S + 
C«f m (t, S )= j Cg B (M)+ for t<8~ (48) 

[ (9nfm)s for t = S 

The two-time correlation function of corresponding global quantities are defined by 

C 9f (t,s):=±J2 C9 J™(t,s), (49) 



(50) 



and the equations of motion for t > s + are given by the following familiar set 

C°f(t,s) + \( 25-1 -S\( C'(t,s)+ 
C«f(t,s) + J [ 35 -3 J [ C*f(t,s) + 

The general solution for these correlation functions in the regime t > s + is given by 

C af (t, s) = A af (s)e- x - (t+s) + E af (s)e- {XM+s) 

+ F af (s)e- {x+t+x - s) + B af (s)e- x+{t+s \ (51) 

C qf (t,s) = \a. (A af (s)e- x - {t+s) + E af (s)e- {x - t+x+s) ) 
5 

+ \a + (F af (s)e^ x+t+x - s) + B af (s)e~ x+{t+s) ) , (52) 
5 

where A a * (s), B a f(s), E u ^ (s), and F a f(s) are arbitrary functions of s. These functions 
are fixed if we require that at t = s the two-time correlation functions become the 
equal-time correlation functions as given in equations (142j) . (145j) . and ( 14*61) . In fact, all 
the arbitrary functions turn out to be independent of s, and are explicitly given by 

A aa = B_ , B a ° = B+ , A aq = -a_5_ , B aq = \a+B+ , 

5 5 

E aa = F aa = ^ ? = J_ a+£ > o ? ^ = l tt _B . (53) 

For these values, the expression for C aq (t, s)+ in equation (151"]) coincides with that for 
C" ?f7 (s,t) + in equation (|5"2"j) . and hence the piecewise solution in t > s extends to the 
t < s~ region, too. 

The equations of motion (l5"0l for the two-time correlations are exact for periodic 
boundary conditions. On the other hand, the system of equations (I4TI) that describes 
the equal-time correlations only holds true in the infinite-size limit L — > oo only when 
at least one the global correlators is much larger than of order 0(1/ L). Therefore, the 
global correlations with the initial state are given by 

C*f(t, 0) = ^ (a'J e" A -' - A a Je~ x+t ^j (54) 
C qf (t, 0) = ^ (a_A*V A -< - a+ Al f e- x ^) (55) 
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where AZ* = a±C af {0) - 5C qf (0). 

We now discuss the behaviour of the global two-time correlators for the same initial 
states as discussed above for the response function. 

(i) For the fully ordered initial states one obtains equilibrium relaxation at the static 
critical point 5=1. We have B- = Bq = 1 and B + = and find 

C^(t, s ) = l + i( e - 2s + e" 2 *) (56) 

such that one has an exponentially fast relaxation towards the global equilibrium 
correlator CZ° = 1. 

(ii) For a fully disordered state the equations of motion do not close and (PUT) are not 
valid. However, the global correlator with the initial state may be read off from 
(I54p and we find C acr (t,0) = 1/L. Similar results may be found for any initial 
temperature T ini > 0. 

(iii) For the partially ordered initial states ■ ■ ■ TTITTITTI ■ " " an d • ■ ■ IITIITIIT • • • we 
find at the critical point 5 = 1 that B- = 1, Bq = —1/3 and B + = 4/9. Hence 

C™(t, s) = 1 - I (e~ 2s + e- 2t ) + - e ~ 2 ^ (57) 
6 9 

which, again, relaxes exponentially towards equilibrium. For the global spin-spin 
correlator with the initial state we find from (154"]) 

C™{t^)= l -- 2 -e- 2t (58) 
which is distinct from the formal s — > limit of the two-time result (E 



For an interpretation in terms of the scaling form (JTJ), recall that C acr (t,s) actually is 
the Fourier transform 

C aa {t, s)= fdr e- iqr C(t, s; r, 0) = s~ b+1/z ( - - 1 V f du F c (u, (59) 

and similarly, a correlation with respect to the initial state is interpreted via (J2J) as 

/ dr e- iqr C(t, 0; r, 0) = t {1 ~ Xc)/z [ du $ c (w) (60) 

JR q =0 JR 

Comparing with the explicit result (158]) for large times, we read off 

Ao = l. (61) 

If we could formally compare eqs. (I57|59l) . we would find b = 1/z = 1/4, but since 
(157]) apparently describes a relaxation towards equilibrium, it is not clear whether the 
non-equilibrium scaling form (159]) is applicable in the present context. 



C aa (t,0) 
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4. Conclusions 

The exact study of one-dimensional kinetic Ising models may provide useful insight into 
the slow relaxation behaviour of many-body systems with strongly interacting degrees 
of freedom[jj] When considering the local spin- flip dynamics given by the rates (J4]), the 
KDH line 7 = 25 is an interesting variant of the celebrated Glauber dynamics 8 — pQ. 
We have generalised the known exact results for the global magnetisation and the global 
three-spin magnetisation [2, H] and have shown how certain global two-time correlators 
and certain global response functions (susceptibilities) may be found, when suitable 
initial conditions are chosen. 

Although the KDH dynamics does satisfy detailed balance, the dynamical behaviour 
close to the critical point at zero temperature is different from the one found for Glauber 
dynamics, since at T = the number of stationary states grows exponentially with the 
number L of lattice sites. This property might be seen as an analogy with the many 
(meta-)stable states in glassy or kinetically constrained systems, such as the Frederikson- 
Andersen model, see e.g. [171 ITS] , or in frustrated magnets, see e.g. [191 EH [21]. In 
particular, analysis of the critical slowing-down near criticality gives the KDH dynamical 
exponent [2] 



in contrast to the result z = 2 of Glauber dynamics [I]. We studied the non- 
equilibrium relaxation of KDH dynamics through the global two-time correlators and 
responses. Although our results appear to be compatible with the usually admitted 
scaling behaviour (T5]) and we have in this way identified some non-equilibrium exponents, 



further tests of non-equilibrium dynamical scaling in the KDH-Ising model are desirable. 
The ageing exponent a could not be determined. Before accepting the formal conjecture 
b = 1/z = 1/4, the applicability of the scaling form ([1]) needs to be checked, which 
cannot be done from the present analytical results alone. These values of exponents 
should be broadly independent of the precise form of the initial states, in agreement 
with the expected universality. We observe that the values of Ac = Xr agree with what 
is found for Glauber dynamics, while the conjectured value of b is different from the 
Glauber dynamics result 6 = 0, see [HI [TDl [TT] . 

Since there is a duality mapping onto a diffusion-annihilation process, it would be 
interesting to see how the results obtained here might bear on that system. In this 
context, we remark that in a similar way as analysed in this paper, one can derive 
closed systems of equations of motion for staggered quantities like M s and T s , but along 
a different line 7 = —25. This line can also be obtained from the line 7 = 25 via a 
gauge transformation o n — > (— l) n a n and 7 — ► —7, which is usually made use to relate 

|| In order to design new slowly relaxing nanosystems, with a view of possible applications to 
information storage, the slow relaxation dynamics in real systems such as the single-chain magnet 
[Mn2(saltmen)2Ni(pao)2(py)2](ClC>4)2 has been explicitly compared to ID Glauber dynamics [16] . 



z = 4 



(62) 



Ac* — — 1? 



(63) 
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ferromagnetic and anti-ferromagnetic Ising models. Graphically, this can be illustrated 
in figure [21 where the 'anti-ferromagnetic' case (—7, 5) can be obtained by a reflection 
about the 7 = line. Quantities such as M s (t) have an immediate interpretation as the 
total particle number in the dual diffusion-annihilation process. 



Appendix: Equations of motion for global correlation functions 

We derive the equations of motion for the correlation functions C^(t) and C g f(t) along 
KDH line 7 = 25. Here, we need not assume that C^(t) = although for r = 1, 2 

it can be easily shown that C[ q (t) = C qa {t) and C£ q (t) = C^it). 

The equation of motion for C^ a (t), for r 7^ 0, is obtained from equation (1151) and 
reads 

^CT(t) = -2C?(t) + 25 (C°Z x {t) + C™ +1 {t)) - $ (C? q (t) + Crit)) . (Al) 

For r = 0, we have C^it) = 1. The dynamic equation for C^ q (t), for r 7^ 0, ±1, 
obtained from equation f|T6|) . is given by 

AcT(t) = - 4C7(t) + 25 (CZM + C^it)) - 5 (CT(*) + c q \t)) 

+ 5 (C^(t) + C^t)) + 6C? A (t) } (A2) 
where C° A (t) := J2 n (a n A n+r ) / L. This is supplemented by 

cro = c™{t) , cg(i) = crit). (as) 

The equation of motion for C qq (t), for r 7^ 0,±1,±2, is obtained from equation (fTTl) 
and is given by 

^CW(f) = - 6C qq (t) + 2*(C£ 1 (f) + C^f)) + 2*(C£i(t) + C&(f)) 

- 5 (C qcr (t) + C?(t)) + 5 (C qA (t) + C Aq (t)) , (A4) 
while that for C% q (t) is obtained from equation (|18p and reads 

—C qq (t) = - 4C qq (t) + 45Cr(t) + 25Crit) - 25C° q it) 
at 

+ 5(c qA (t) + C Aq (t)y (A5) 

For r = 0, 1, we have 

C qq (t) = l , C qq {t)=Cr{t). (A6) 

The equations of motion for C qa , C qq } are not closed as they also involve C Aq 

terms. Hence explicit solutions can only be found via approximate methods, such as 
mean-field theories or truncations. 

However, one may obtain the equations of motion for the global correlation functions 
C 9 ^{t). In the expansion 



d 



^^EE^M- (AT) 



dt w L 2 ^ dt 
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upon substituting equation f)15p and rewriting in terms of correlation functions, we 
obtain 

^-C aa (t) = (AS - 2)C aa (t) - 25C qa (t) + j(l-25C™{t)+5C™{t)). (A8) 
Similarly, in the expansion 

s c "w = fE E ^^> + i(s c «"w + 2 ^'w) < A9 > 

upon using equations (|T6l) . (1A3I) . and (1A1[) . we find 

Ac^(t) = 35C" T<T (t) + (25 - 4)C CT9 (t) - SC qq (t) 
at 

+ + 2(1-5)CTW + (l-25)C™(t) + 26C 3 TCT (t) - 25C° 9 (*)). (A10) 
Finally, from the expansion 

s c "«^E E s^ + Ks^w+s^w) (A11) 

and upon using equations (TT7T) . ( 1A6I) . ( 1A5j) . and ( 1A1I) . we obtain 

^C"(t) = 65C CT9 (t) - 6C"(t) (A12) 

+ v(3 - 25Cr(0 - ^2 CT (t) + 4C 3 ff<T (t) + 25CT(*) 
- 46C?{t) - 26(C° q (t) + Cf (t)) + 2C*(f)) • 

The equations of motion for {C™, C"" 7 , C qq } also involve {Cf, C 2 CTCT , C™, Cf, C^ 9 , 
Although these equations are not closed in general, they sometimes do close in the ther- 
modynamic limit, at least for certain initial conditions. 

For example, in the case of a fully disordered initial state, C 9 ' are of order 
0(l/y/L), while C gf are of order 0(1/ L). Now the solutions Cf(t) will be of the 
form (t) + Q g r f (t)/VZ. This set of functions {f 9 r f (t)} is needed to compute C 9f (t) to 
order 0(1/ L), but for which no closed system of equations is known. On the other hand, 
if we choose an initial state such that C 9 -f(0) = 0(1), the limit L — > 00 in equations ( IA81 
IA10t IA12I) can be taken and the resulting closed system fj4"Tj) may be solved explicitly, 
as discussed in the main text. 
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